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Stochastic  Differential  Equations  for  Power  Law  Behaviors 

Bo  Jiang  Roger  Brockett  Weibo  Gong  Don  Towsley 


Abstract — We  present  in  this  paper  some  simple  stochastic 
differential  equations  that  leads  to  lower-tail  and/or  upper 
tail  power  law  behaviors.  We  also  present  a  model  with  bi¬ 
directional  Poisson  counters  that  exhibits  power  law  behavior 
near  a  critical  point,  which  might  be  of  interest  to  statistical 
physics. 

I.  Introduction 

As  has  been  observed  many  times  in  recent  years,  empir¬ 
ical  studies  dealing  with  a  variety  of  subject  matter  produce 
data  showing  power  law  histograms  extending  over  several 
decades  or  more.  The  diversity  of  application  areas,  and 
the  extent  of  the  data  have  inspired  researchers  to  look  for 
general  principles  that  would  make  it  possible  to  trace  the 
individual  phenomena  back  to  one  or  more  common  features. 
Among  recent  work  along  these  lines  we  call  attention  to  [?], 
[?],  [?],  [?].  In  this  paper  we  pursue  the  idea  that  certain 
simple  forms  of  first  order  stochastic  differential  equations 
have  steady  state  densities  which  show  lower  tail  and/or 
upper  tail  power  law  behaviors,  depending  on  the  values  of 
the  parameters.  To  the  extent  that  the  form  and  parameters 
of  a  differential  equation  are  often  more  readily  identified  in 
the  modeling  process,  the  differential  equations  we  describe 
may  be  thought  of  as  providing  a  more  direct  explanation  of 
power  law  behavior. 

In  their  interesting  paper  [?],  Reed  and  Hughes  consider 
the  probability  density  of  the  value  of  a  growing  exponential 
sampled  at  an  exponentially  distributed  random  time.  As  they 
observe,  this  is  easily  seen  to  give  a  power  law  distribu¬ 
tion.  Here  we  consider  a  different  situation  involving  the 
steady  state  density  associated  with  a  stochastic  differential 
equation.  The  equation  describes  a  situation  in  which  the 
quantity  of  interest  decays  to  zero  following  an  exponential, 
X  =  —aX,  but  is  incremented  by  a  fixed  amount  a  at 
random  times,  the  times  having  an  exponential  distribution. 
We  show  that  for  a  range  of  parameter  values  the  steady  state 
distribution  of  A'  exhibits  a  power  law  lower  tail.  The  fact 
that  we  deal  with  the  steady  state  property  of  an  ongoing 
dynamics  gives  our  work  a  different  set  of  possibilities  for 
interpretation.  The  basic  reason  as  to  why  the  distribution  of 
values  in  steady  state  has  the  same  form  as  the  distribution 
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of  values  obtained  by  sampling  at  a  random  time,  lies  in 
the  form  of  the  drift  term  associated  with  the  differential 
equation. 

When  considering  lower  tail  behavior,  we  find  it  useful  to 
distinguish  between  what  might  be  called  “two  sided”  lower 
tails  and  “one  sided"  lower  tails.  The  latter  may  occur  when 
dealing  with  an  intrinsically  nonnegative  quantity,  such  as 
queue  length  whereas  the  former  might  apply  in  a  situation  in 
which  the  histogram  extends  in  both  directions  from  a  critical 
value,  as  might  be  the  situation  for  some  types  of  populations 
near  a  phase  transition.  Our  models  for  one  sided  lower  tail 
behavior  focuses  on  the  steady  state  behavior  associated  with 
stochastic  differential  equations  containing  a  Poisson  counter 
N  and  taking  the  form 

dA  =  f(X)dt  +  g{X)dN. 

For  two  sided  lower  tails  we  consider  slightly  more  compli¬ 
cated  models  of  the  form 

dA  =  f(X)dt  +  gi(X)dNi  -  g2(X)dN2. 

We  also  study  the  more  popular  case  of  upper  tail  power 
law.  A  simple  transformation  is  shown  to  convert  one  sided 
lower  tail  power  laws  into  upper  tail  power  laws.  In  view 
of  dynamic  systems,  evolutions  of  this  type  can  also  give 
generative  interpretations  of  many  power  law  upper  tails 
observed  in  real  data. 

In  Section  ??  we  first  present  our  motivating  example 
of  a  simple  Poisson  counter  driven  SDE,  the  steady-state 
distribution  of  which  was  shown  by  Brockett  [?]  to  exhibit 
power  law  behavior  near  the  origin.  We  then  introduce  a 
simple  transformation  to  develop  a  similar  SDE  that  leads  to 
an  upper  tail  power  law.  We  observe  that  various  distributions 
can  be  generated  via  simple  modifications  of  the  “drift 
term”  in  the  SDE.  In  Section  ??,  we  add  a  Brownian 
motion  term  to  obtain  a  similar  result  as  Reed  [?].  All  of 
these  cases  demonstrate  that  “random  multiplication  with 
exponential  stopping  time  will  lead  to  power  law  behaviors”. 
In  Section  ??  we  develop  an  SDE  driven  by  Poisson  counters 
in  both  positive  and  negative  directions.  We  show  that  the 
steady-state  distribution  can  exhibit  power  law  behavior  near 
a  critical  point.  This  may  have  implications  in  statistical 
physics  since  a  discontinuity  occurs  in  a  surprising  way. 
Section  ??  concludes  the  paper. 

II.  SDE  DRIVEN  BY  A  POISSON  COUNTER 

The  following  SDE  was  considered  by  Brockett  [?], 

dAt  =  —aXtdt  +  crd/Vt, 


(1) 


where  a,  a  >  0  and  N  is  a  Poisson  process  of  intensity 
A.  By  Theorem  6.2  in  [?],  there  is  a  unique  adapted  RCLL 
process  {Xt}  satisfying  (??)  and 

sup  E[X2]  <  oo  (2) 

te[o,T] 

for  any  T  £  [0,  oo).  Similar  arguments  for  the  existence 
of  solutions  apply  to  all  the  other  SDE’s  considered  in  this 
paper  and  will  not  be  repeated. 

Now  let  ipk(x)  =  elkx ,  k  €  R.  By  Ito’s  formula. 


d tpk{Xt)  =  -  iakXtipk(Xt)dt 

+  bPk(Xt-  +  cr)  —  if>k(Xt-)]dNt 

=  -  akd4’kdft]  dt  +  (e*k°  -  \)MXt-)dNt. 

Taking  expectation, 

A  A 

^<MM)  =  -ak—$x(k,t)+\(eik°-l)<S>x(k,t),  (3) 

where  &x(k,t)  =  E [tpk(Xt)\  is  the  characteristic  function 
of  Xt  and  the  change  of  the  order  of  differentiation  and  ex¬ 
pectation  can  be  justified  by  (??)  and  Lebesgue’s  Dominated 
Convergence  Theorem.  Equation  (??)  can  be  solved  by  the 
method  of  characteristics  (see  e.g.  [?])  to  yield 


<l>x  ( k,t )  =  ( ke  at ,  0)  exp 

After  a  change  of  variable  u  = 


gikc Tea^s  O  _ 


>0  L 


ds 


(4) 

a  exp(a(s— t)),  (??)  becomes 


{k,  t )  =  <f>x  (ke  at ,  0)  exp 


yielding 


&x(k,  oo)  =  exp 


piku  _  i  > 

- -du\ 

u  J 


This  result  can  also  be  obtained  as  a  consequence  of  Theorem 
2  of  [?].  By  Lemma  53.2  of  [?],  the  steady-state  distri¬ 
bution  is  absolutely  continuous  with  continuous  density  on 
(0,  oo).  Moreover,  the  density  is  continuously  differentiable 

on  (0,  oo)  if  A  >  a. 

Note  that  <5>x(k,  oo)  satisfies  the  following  equation,  ob¬ 
tained  from  (??)  by  setting  the  right-hand  side  to  zero, 

-ak^-$(k)  +  X(eika  -  l)$(jfe)  =  0. 

d  k 

Thus  the  corresponding  density  fx(%)  satisfies 


a—[xf(x)}  +  A  f(x  -a)-  A  f(x)  =  0. 


(6) 


It  was  shown  in  [?]  that  fx(x)  =  0  for  x  <  0.  We  can  arrive 
at  the  same  conclusion  in  a  more  intuitive  way  by  examining 
(??).  Indeed,  note  that,  in  spite  of  its  initial  value,  X  will 
eventually  become  positive  and  remain  so  from  that  point 
on.  Now  using  fx(x)  =  0  for  x  <  0,  (??)  can  be  solved 
recursively  to  give 


fx(x)  =  Cx«  \  x  £  (0,  cr] , 


and,  for  x  £  ( na ,  na  +  a],  n>  1, 


fx(x)  =  fx(ncr) 


\na/ 


■u  fx(u-cr)du, 


where  the  constant  C  is  determined  by  the  normalization 
condition 


/OO 

fx(x)dx  =  1. 

-OO 


Note  that  fx(x)  has  a  power  law  at  its  lower  tail. 

We  briefly  mention  that  (??)  can  have  the  following 
generalization. 


dXt  =  AXidt  +  bdNt, 


where  X  is  an  Revalued  process,  A  an  n  x  n  stable  matrix 
and  b  £  Rn .  Then  the  steady-state  distribution  has  the 
following  characteristic  function 

<I>x(k,  oo)  =  exp  |  A  j  [exp  (ikT  eAsb)  —  l]  ds|  , 

where  k  £  Wn,  and  the  projection  of  X  onto  a  left 
eigenvector  of  A  exhibits  power  law  near  the  origin. 

Now  we  introduce  the  simple  transformation  Yt  =  Xt  1 
to  convert  the  lower  tail  power  law  into  an  upper  tail  power 
law.  For  y  >  £  =  cr-1,  the  steady-state  density  of  Y  is 

fv(y)  =  fx(y^)y~2  =  Cy y  £  [e,oo). 

Using  Ito’s  formula,  we  can  get  the  equation  governing  the 
evolution  of  Yt, 

Yt 2 

dYt=aYtdt - ~^—dNt,  (7) 

£  +  Yt- 

Note  that  at  each  jumping  point,  Y  drops  to  ,  which 

depends  on  Yt  _ . 

We  can  modify  the  coefficient  in  front  of  d Nt  in  (??) 
so  that  it  always  restores  the  process  to  a  fixed  point.  The 
equation  then  becomes 

dZt  =  aZtdt  +  Oo  -  Zt_)dIV(,  (8) 


with  the  corresponding  equation  for  the  characteristic  func¬ 
tion  &z(k,t)  of  Zt  being 

A  A 

~^z(k,t)  =  ak—^z(k,t)  -  X^z(k,t)  +  Xeikzo.  (9) 
at  ok 

Equation  (??)  can  be  solved  again  by  the  method  of  charac¬ 
teristics,  yielding 

$z(M)  =e~xt<$>z(keat,  0) 

+  A  J  exp  |  A(s  —  t)  +  izoke °4t-s)  |  ds. 

(10) 

After  a  change  of  variable  z  =  Zoea^  s\  (??)  becomes 

at  A  i 

A  fz°e  f  z  \  ~  “_1 

*z(k,t)  =  e~xt<f>z(keat,  0)4- -  /  -  eikzdz, 

azo  JZo  \z0J 

from  which  we  can  read  off  the  distribution  function, 

Fz(z,t)  =  e-XtFz(ze-at,0)  +  (1  -  e~xt)G(z,t), 


where  G(z.  t)  is  a  truncated  Pareto  distribution. 


where 


G(z,t) 


As  t  — >  oo,  the  distribution  Fz(z,t ) 
distribution 


2  <  Z0 , 

zo  <  z  <  z0eat, 
z  >  z0eat. 
approaches  a  Pareto 


Fz(z,  oo )  =  1  - 


z  >  z0- 


A  closely  related  model  of  deterministic  exponential 
growth  with  exponential  stopping  time  was  analyzed  in  [?]. 
We  note  by  passing  that  the  proportional  growth  is  critical 
in  generating  power  law.  Had  the  growth  term  in  (??) 
been  a.Zfdt  for  some  S  £  [0, 1),  the  resulting  distribution 
would  have  been  a  left  truncated  Weibull  distribution  with 
distribution  function 

Fz(z,  oo )  =  1-exp  j  — ^  (z1-6  -  j  ,  z>  z0. 


III.  SDE  DRIVEN  BY  BOTH  BROWNIAN  MOTION  AND 
Poisson  Counter 

In  this  section,  we  add  a  Brownian  motion  component  to 
(??),  which  becomes 

dXt  =  ,iXtdt  +  aXtdWt  +  {x0  -  Xt- )d Nt,  (11) 


where  /x,  Xq  £  R,  cr  >  0,  W  is  a  standard  Brownian  motion 
and  N  is  a  Poisson  process  with  density  A,  independent  of 
W.  This  is  a  geometric  Brownian  motion  with  Poisson  jumps 
which  always  reset  the  motion  to  a  fixed  state  Xq .  A  similar 
model  was  analyzed  in  Reed  [?]. 

Let  Yt  =  \ogXt  and  yo  =  log^o-  Then  Ito’s  formula  gives 

d Yt  =  (y  -  2  j  dt  +  adWt  +  ( y0  -  Yt_)dNt, 


which  is  a  Brownian  motion  randomly  reset  to  yo  by  Poisson 
jumps.  Let  V’fe(y)  =  ezfey  as  in  the  previous  section.  By  Ito’s 
formula. 


d  MYt) 


iki/jk(Yt ) 


dt  +  odWt 


-  1 a2k2MYt)dt  +  (eikyo  -  ^k{Yt-))dNt. 


Taking  expectations,  we  get 


8 ' 
dt 


$v(M) 


k - cr2k2  —  A 

2 


$y(M) 


+  \elky\ 


where  <l>y(k.  t)  is  the  characteristic  function  of  Yt.  The 
solution  is 


<l>y(fc,t)  =  $y(fc,Oo) 

+  e-At[<Ly(fc,0)  -  <S>Y{k,oo)]el{yt--*a2t)k-^2tk2 , 


<I>y  (fc,  oo)  = 


—\eikyo 


i(fi  —  1 a2)k  —  \cr2k2  —  A 


Now  we  can  find  the  steady-state  density  of  Yt  as  t  — >•  oo 
by  taking  the  inverse  Fourier  transform  of  <by(fc,  oo), 


f  jYLeP(y-yo) 

fy(y)  =  aJl  {y_y’o) 

{ct+pe  ’ 


y  <  yo, 
y  >  yo, 


where  a  >  0  and  —  (3  <  0  are  the  two  roots  of  the  following 
quadratic  equation. 


1 


+  I  /x  -  I  7  -  A  =  0. 


1 


(12) 


Going  back  to  X,  we  get  the  steady-state  density  of  Xt  as 

t  — >  oo 


fx{x)  =  fY(logx)x  1 

\  x0  >  cce(0,a;o], 

=  0  a+aP)X0(-aJ1  [  J  (13) 

1*0^  (^)  >  x£[x0,oo). 

which  is  the  double  Pareto  distribution  of  Reed  [?]. 

Motivated  by  the  connection  between  (??)  and  (??),  we 
also  consider  the  following  SDE 

Z2 

dZt  =  pZtdt  +  aZtdWt  -  *-  dJVt,  (14) 

Zt-  +£ 

with  Z0  >  0.  Let  Ut  =  Z^ 1  •  Then,  by  Ito’s  formula, 

dUt  =  -(/x  -  a2)Utdt  -  aUtdWt  +  e_1d Nt. 

By  the  same  procedure  as  before,  we  get  the  equation  for 
the  characteristic  function  &u(k,t)  of  U , 

o  o 

—  <5>u(k,t)  =  -  {n  -  a2)k—<Pu{k,t) 

+  ^°'2fc2(|^2  $u(M)  +  (e*fc/£  - 

which  is  the  Fourier  transform  with  respect  to  the  variable 
y  of  the  following  Fokker-Planck  equation  for  the  density 
of  Yt, 

^ fv(y,t )  =  O-  °2)-j^[yfY{y,t)\  +  *o-2 ^[y2 fY(y,t)} 
+  A/y(t/-e_1,i)  -  XfY{y,t). 

In  steady  state,  the  density  fY(y)  satisfies 

(j*  ~  ^-^[yfriy)]  +  b2/v(y)] 

+  A fY{y  -  £_1)  -  A fY(y)  =  0. 

For  y  £  (0,£_1],  this  reduces  to 

(a*  — cr2)^ l-[yfY(v)\  +  \a2^X2  iy2fy(y)\  -  =  °- 

(15) 

The  general  solution  to  (??)  is  given  by 

fY(y)  =  Cy*-1  +  Dy-13-1,  ye(0,e_1], 


where  a  and  j3  are  as  in  (??).  The  integrability  of  fy  (jj) 
requires  that  D  =  0,  so 

fy{y)  =  Cya~1,  ?/e(o,e_1]. 

Therefore, 

fx(x)  =  fY(x^1)x~2  =  Ca;_a_1,  x  £  [e,oo), 

which  has  the  same  upper  tail  power  law  exponent  as  in  (??). 
This  is  intuitive  since  the  difference  of  the  two  models  lies 
in  the  range  of  small  x.  In  fact,  both  (??)  and  (??)  have  the 
same  general  form 

dXf  =  gXtdt  +  aXtdWt  +  \g(Xt_)  -  Xt_]d Nt,  (16) 

where  0  <  g(-)  <  B  is  bounded.  The  density  of  Xt  satisfies 
the  same  Fokker-Planck  equation  for  x  £  (B,  oo), 

= ~L{ijxfx) + l^{a2x2fx)  - xfx- 

Therefore,  the  upper  tail  of  the  steady  state  density  has  the 
same  functional  form,  independent  of  g(-).  By  taking  the 
limit  er  — >  0,  we  observe  that  the  remark  here  also  applies 
to  (??)  and  (??). 

We  note  also  the  possibility  of  extending  (??)  to  the 
multivariate  case, 

dXi=IMXidt+aiXidWi  +  YJ[9ijOL-)-Xi-}dNj,  (17) 

3 

where  we  have  omitted  the  subscript  t.  As  in  (??),  each 
Xi  has  a  power  law  upper  tail.  Thus  (??)  generates  a  set 
of  correlated  random  variables  with  power  law  marginal 
distributions. 

IV.  SDE  DRIVEN  BY  BI-DIRECTIONAL  POISSON 
COUNTERS 

In  this  section,  we  consider  the  following  SDE  driven  by 
bi-directional  Poisson  counters. 


This  result  can  also  be  obtained  using  Theorem  17.5  of  [?]. 
Lemma  2  of  [?]  shows  that  &x(k,  oo)  is  the  characteristic 
function  belonging  to  an  absolutely  continuous  distribution 
and  the  density  is  continuous  if  and  only  if  Ai  +  A2  >  a.  If 
Ai  =  A2  =  A  and  or  =  172  =  cr,  the  steady-state  distribution 
will  be  symmetric  around  the  origin.  If,  in  addition,  a  = 
<To  A- 2,  then  as  A  — >  oc,  (??)  converges  to  the  Ornstein- 
Uhlenbeck  process  and 

T>x(fc,  00)  -A  exp  . 

i.e.  the  characteristic  function  of  the  normal  distribution 

2 

A/”(0,  ^),  as  expected. 

Setting  the  right-hand  side  of  (??)  to  zero,  we  get  the 
differential  equation  satisfied  by  $x{k,  00), 

ak^x(k,t)  =  [A1(e*fecri  -  1)  +  A2(e-ifeTa  -  1)]  <MM)- 
Thus  the  steady-state  density  satisfies  the  following  equation, 

a^-[xf(x)]+X1f(x-a1)-(X1+X2)f{x)+X2f{x+a2)  =  0. 

(20) 

Figure  ??  shows  some  steady-state  densities  for  (??) 
obtained  from  simulation,  where  we  have  set  or  =  a2  =  1 
and  Ai  =  A2.  Note  that  as  Al^Aa  becomes  larger,  the  density 
becomes  smoother,  consistent  with  Lemma  2  of  [?].  As 
Al+A2  becomes  smaller,  the  density  becomes  more  sharply 
concentrated  around  zero.  In  the  case  A1  +  A2  <C  a,  f(x—a\) 
and  f(x  +  a 2)  are  negligible  compared  to  f(x)  for  small 
x  and  hence  (??)  can  be  approximated  by  the  following 
equation, 

a~[:c/(:r)]  -  (Ai  +  A 2)f(x)  =  0. 

which  can  then  be  solved  to  give 


dXt  =  a(g  -  Xt)dt  +  <JidPt  -  a2dNt, 

where  /i  £  R,  a,<Ji,cr2  >  0,  and  P,  N  are  two  independent 
Poisson  processes  with  intensities  Ai,  A2,  respectively.  By  a 
simple  shift  of  the  origin,  we  may  assume  without  loss  of 
generality  that  g,  =  0  and  the  equation  then  becomes 


dXt  =  —aXtdt  +  <J\dPt  —  a2dNt.  (18) 

The  characteristic  function  (\>x  (k,  t)  satisfies  the  following 
equation, 

o  a 

—  $x(k,t)  =-ak—$x{k,t) 

+  [A1(eife<Tl  -  1)  +  \2(e~ika2  -  1)]  <MM), 

(19) 

the  solution  to  which  is 

$x(M) 


-at 


=<&x{ke  ,  0)e  “  J<Ti 

Therefore, 


t=idu-^  f~pe  Qt  ?Xix=±du 


-(k,  00)  =  e 


Ar  T?  1  2  ro  <XXzldu 

-  S)  CX  JU  1L  OC  J  — (TO  'LL 


f(x)  =  C\x\  2  1,  0  <  \x\  <C  1. 


Figure  ??  plots  the  steady-state  densities  in  log-log  scale 
with  the  reference  lines  of  slope  Al+A2  —  1  superimposed. 
The  approximation  is  very  good  near  the  origin. 

This  approximation  can  be  made  more  rigorous  using 
Theorem  53.8  of  [?].  We  will  present  a  simplified  and  more 
straightforward  analysis  of  the  behavior  of  the  density  near 
the  origin  when  Ai  +  A2  <  a.  Let 


and 


$1  (k)  =  exp 


„iku  _  1  'v 

- - -du[ 

u  J 


$2  (k)  =  exp 


piku  __  1  > 

- - -du  \  , 

u  J 


which  are  the  characteristic  functions  of  two  absolutely 
continuous  distributions.  Denote  their  densities  by  g(x)  and 
h(x),  respectively.  Then  [?]  shows  that  g(x)  has  support 
on  [0,  cxj)  and  g(x)  =  CxXr~Y  for  x  £  (0,  cti].  Similarly, 
h(x)  has  support  on  ( — 00, 0]  and  h(x)  —  D \x\yy~1  for 


Let  m  =  minjfT!,  a2j  and  e  €  (0,  m).  For  0  <  x  <  m  —  e, 
we  have 

pOO 

f(x)  =  /  g{y)h(x  -  y)dy 

J  X 

f-m  x  x  roo 

=  /  Cy^^1D(y  ~  x)^^1dy  +  /  g(y)h(x  -  y)dy 
J  x  J  m 

=  CDxX-^-1  j™  '(«  •  1)*  'du 

pOC 

+  /  g(y)h(x-y)dy 
J  m 

Aj+A2  -j 

=  A\X  a  +  A2 , 

where  A\  =  CD  Jp  u^-1(w  —  l)~^'~1du  and  A2  = 
fm  9{y)h{x—y)dy.  It  is  shown  in  [?]  that  g(x)  is  continuous 
on  R.  \  {0}  and  g'(x)  <  0,  so  it  is  uniformly  bounded  on 
[m,  00)  and  hence 

pOO 

A2  <  sup  g(y)  /  h(x  -  y)dy  <  sup  g(y)  <  00. 
y>m  J  m  y>m 

If  Ai  +  A2  <  a,  then 


A 1  >  CD  /  u  «  1(u  —  1)  “  1du  >  0 


and 


/°°  \  \ 

—  l)^"1du  <  00. 

A  similar  analysis  applies  for  the  case  x  £  (— m  +  e,0). 
Therefore, 

f(x)  =  0  ,  as  x  — >  0. 

If  Ai  +  A2  =  cr, 

Ai  =  CD  (  *  u ^(u  -  1)  «  'du 


>  CD  /  u  “  u  “  du  =  CD  log  — , 
J 1  x 

and  for  x  < 

.4,  =  CD  jj  ««  '(«  1)*’!  'du 

<  CD  f  (u-  l^^du  +  CD  f*  tM-1 

=  CD^-  +  CD21~xx  log 
A2 

A  similar  analysis  applies  for  x  <  0.  Therefore, 
f(x)  =  0(—  log  |x|),  as  x  — >  0. 


Xl-l 


du 


V.  Conclusions 

We  presented  some  simple  stochastic  differential  equations 
that  lead  to  lower  tail  and/or  upper  tail  power  law  behaviors. 
Some  of  the  results  are  known  but  the  derivations  are 
different.  We  also  presented  a  model  with  two  opposite 
Poisson  counters  and  an  exponential  decaying  term.  This 
model  exhibits  power  law  behavior  near  a  critical  point, 
which  might  be  of  interest  to  statistical  physics. 


